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Abstract 

o 

^ ■ Here, cosmology of the late and future universe is obtained from 

f(R)— gravity with non-linear curvature terms R 2 and R 3 (R being 
the Ricci scalar curvature). It is different from f(R)— dark energy 
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models [6] where non-linear curvature terms are taken as gravitational 
alternative of dark energy. In the present model, neither linear nor 
non-linear curvature terms are taken as dark energy. Rather, dark 
energy terms are induced by curvature terms in the Friedmann equa- 
tion derived from f(R)— gravitational equations. It has advantage over 
f(R)— dark energy models in the sense that the present model satisfies 
WMAP results and expands as ~ t 2 ^ 3 during matter dominance. So, 
it does not have problems due to which f(R)— dark energy models are 
criticised in [7] . Curvature- induced dark energy, obtained here, mimics 
phantom. Different phases of this model, including acceleration and 

deceleration during phantom phase, are investigated here. It is found 
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that expansion of the universe will stop at the age (3.87to + 694.4)kyr 
(to being the present age of the universe) and, after this epoch, it will 
contract and collapse by the time (336.87t + 694.4)kyr. Further, it is 
shown that universe will escape predicted collapse (obtained using clas- 
sical mechanics) on making quantum gravity corrections relevant near 
collapse time due to extremely high energy density and large curvature 
analogous to the state of very early universe. Interestingly, cosmologi- 
cal constant is also induced here, which is extremely small in classical 
domain, but very high in the quantum domain. 

Key Words : f(R)— gravity cosmology, phantom dark energy, ac- 
celeration and deceleration in phantom universe, contraction of the 
universe, cosmic collapse in future and quantum escape of collapse. 

1. Introduction 

Astrophysical observations, made at the turn of the last century [1, 2] 
show conclusive evidence for acceleration in the late universe, which is 
still a challenge for cosmologists. Theoretically, it is found that dark 
energy (DE) violating strong energy condition (SEC) or weak energy 
condition (WEC) is responsible for it. So, in the recent past, many DE 
models violating one of these conditions, were proposed to explain the 
late cosmic acceleration taking some exotic matter being the possible 
DE source. A comprehensive review of these models is available in 
[3] . As the exotic matter is not known, different scalar fields with spe- 
cial attention to quintessence, phantom, K-essence and tachyons, were 
tried upon to probe actual nature of the matter resposible for required 
DE [3] . Apart from, field-theoretic models, some fluid dynamical mod- 
els were also proposed. Chaplygin gas and generalized Chaplygin gas 
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models of DE are prominant among these, as it has super-symmetric 
generalization and negative pressure needed for violation of SEC or 
WEC[3, 4]. 

As no satisfactory model came up, it was thought that curvature, 
being the source of gravitation, could also be a source of DE. The 
Einstien-Hilbert term R/l67rG (R being the scalar curvature and G 
being the gravitational constant) yields the geometrical component in 
Einstein's gravitational equation, so it was realized to use non-linear 
terms of curvature replacing matter lagrangian in the Einstein's grav- 
itational action. Motivated by this idea, S. Capozziello et al. and 
S.M. Carroll et al. proposed non-linear terms of curvature R~ n with 
n > as a possible source of DE [5]. Although this model explained 
late cosmic acceleration, it exhibited instability and failed to satisfy 
solar system constraints. Further, this idea was taken up by Nojiri and 
Odintsov and models for gravitaional alternative of DE were proposed 
taking different forms of f(R) other than R~ n . These improved models 
satisfied solar system constraints exhibiting late cosmic acceleration for 
small curvature and early inflation for large curvature. Thus, in f(R) — 
dark energy models, non-linear curvature terms were considered as an 
alternative for DE [6, for detailed review]. Recently, Amendola et al. 
have criticized f(R)— dark energy models and have shown that these 
models with dominating powers of R for large or small R can not yield 
viable cosmology as results contradict the standard model and do not 
satisfy Wilkinson Microwave Anisotropy Probe (WMAP) observations, 
though these models pass solar system constraints and explain late ac- 
celeration [7]. Moreover, it is also shown that the most popular model 
with f(R) = R+ aR m + f3R~ n (m > 0,n > 0) considered in [6] is unable 
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to produce matter in the late universe prior to the beginning of late 
acceleration [7]. 

The recent developements against f(R)— dark energy models indi- 
cate the need to think for some alternative approach to get DE from 
curvature. Such an approach has been proposed in [8, 9, 10, 11, 12] 
and it is adapted here too. This apporach yields a cosmology emerging 
from f(R)— gravity satisfying WMAP results as well as results of the 
standard cosmological model. On the contrary, these results (WMAP 
results as well as results of the standard cosmological model) are not 
staisfied by f(R)— dark energy models well as it is noted in [7]. There is 
a crucial difference between the cosmology from f(R)— gravity, being 
addressed here as well as in [8, 9, 10, 11, 12] and f(R)— dark energy 
models. In the latter case, non-linear terms of curvature are taken as 
dark energy lagrangian a priori. On the contrary, in the former case, 
neither linear nor non-linear term is considered as dark energy. Rather, 
in the present cosmological model from f(R)— gravity, DE is induced 
by linear (Einstein-Hilbert term) as well as non-linear terms R 2 and 
R( 2 + r ) in the action ( r being a real number such that (2 + r) > 0). In 
f(R)— dark energy models, dark energy terms depend on f(R) terms 
and its derivative F = df/dR. In the former case, induced DE terms 
depend on the scale factor a(t) of the homogeneous and flat of Fried- 
mann - Robertson - Walker (FRW) model of the universe. 

In what follows, f(R) — gravity with R 2 and R 3 is considered. It 
is found that curvature terms induce dark energy, dark matter and 
cosmological constant in the Friedmann equation (FE) for the late uni- 
verse, derived from f(R)— gravitational equations. It is interesting to 
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see that curvature-induced dark energy, obtained here, mimics phan- 
tom with equation of state (EOS) parameter w = —5/4. Moreover, 
FE contains phantom DE term as Pde[l — Pde/2A]. The correction term 
—p\ e /2\ with A being the cosmic tension [10, 12] is analogous to such 
term in RS-II model based FE [13] as well as loop quantum gravity 
correction [14]. Like references [10, 12], here also, this term is obtained 
from f(R)— gravity. Here, cosmic tension A is evaluated to be 5.77p° r 
with p° r being the present critical density of the universe. Further, it 
is shown that universe, derived by curvature-induced dark matter, de- 
celerated upto time 0.59t - At this epoch and small red-shift z = 0.36, 
transition from deceleration to acceleration took place. Interestingly, 
it is noted that as phantom energy density increases, effect of the term 
—p\ e /2\ gradually increases. As a result, it is found that universe will 
super-accelerate (expansion with high acceleration) during the period 
0, 59t < t < 2.42t , it will accelerate (expansion with low accelera- 
tion) during the period 2.42to < t < 3.44to and, during the period 
3.44t < t < 3.87to, universe will decelerate even during the phantom 
phase. Phantom-dominance will end when p = 2A = 11.54p° r at time 
t = 3.87to causing re-dominance of dark matter with the re-beginning 
of decelerated expansion. It is natural to think that, even during this 
phase with deceleration (derived by dark matter), phantom DE will 
grow with expansion causing pde[l _ Pde/2A] negative due to pdc > 2A 
during re-dominabce of dark matter. In the beginning of this phase, 
this term can be ignored, but its effect will grow causing expansion 
to stop at time t m = 3.87t + 694.4kyr. As a consequence, universe 
will retrace back and contract. Further, it is found that by the time 
i co \ = 336.87to + 694.4kyr, universe will collapse and energy density will 
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diverge. It is argued that this state of the future universe is analogous 
to the state of very early universe, where energy density is extremely 
large and curvature is very high. So, like early universe, it is reason- 
able to think that quantum gravity effects will be dominant near the 
collapse time. Motivated by this idea, quantum corrections are made 
in FE in the small interval of time around the time of collapse and it 
is found that cosmic collapse, predicted by classical mechanics, can be 
avoided. Earlier too, Quantum gravity corrections were done in the 
future phantom universe to avoid finite time singularities [15, 16, 17]. 
The paper is organized as follows. 

In section 2, phantom-phase of the late and future universe is in- 
vestigated. Transition from deceleration to acceleration is probed and 
cosmic tension is evaluated. Section 3 contains discussion on energy 
conditions during phantom era and its cosequences. In section 4, re- 
dominance of dark matter, maximum expansion, contraction and pos- 
sible cosmic collapse using the classical approach have been addressed. 
Avoidance of collapse, using quantum gravity effects, is demonstrated 
in sub-section 4(c). Important results are summarized in the last sec- 
tion. 

Natural units(A;£ = h = c = 1) (where ks, h, c have their usual 
meaning) are used here. GeV is used as a fundamental unit and we 
have lGeV" 1 = 6.58 x lCT 25 sec. 

2. Phantom phase of the late universe and future universe 



from f(R) -gravity 



Here, the action for f(R)— gravity is taken as 




(2.1) 
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where G = Mp 2 [Mp = 10 19 GeV is the Planck mass), a is a dimension- 
less coupling constant and (3 is a constant having dimension (mass) -2 
(as R has mass dimension 2). 
Action (1) yields field equations 

Y^q( R ^ - \d^ R ) + a(2V^V,i? - 2g^DR - -g^R 2 + 2RR^) 

+[3(3V^ U R 2 - 3g^nR 2 - ^g^R 3 + 3R 2 R^) = (2.2a) 
using the condition SS/Sg^ = 0. The operator □ is defined as 

Taking trace of (2.2a), it is obtained that 

-— ^- - 6(a + 3f3R)UR - 18/3Vi2V„i2 + f3R 3 = (2.3) 
167rG 

In (2.3), (a + 3/3 R) emerges as a coefficient of OR due to presence 
of terms aR 2 and PR 3 in the action (2.1). If a = 0, effect of R 2 
vanishes and effect of R 3 is switched off for /3 — 0. So, an effective 
scalar curvature R is defined as 

-fR = a + 3pR, (2.4) 

where 7 is a constant having dimension (mass)" 2 being used for dimen- 
sional correction. 

Connecting (2.3) and (2.4), it is obtained that 

^ 1 uf> s 1 r 1 « 2 i 1 r 1 « 2 i 1 

I 1 d V J?V J? 

R 11 "67 116^67 3/?J 67Ll67rG 9/?J Tj R 

R 

■foR-3a\. (2.5) 



54/3 1 

Experimental evidences [18] support spatially homogeneous flat model 
of the universe 

dS 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 } (2.6) 
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with a(t) being the scale- factor. 

In this space-time, (2.5) is obtained as 

~ a ~ 
-R-3-R 
a 
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[yR - 3a] 



(2.7) 



due to spatial homogeneity. 

For a(t), being the power-law function of t, R ~ a~ n . For example, 
R ~ a~ 3 for matter-dominated model. So, there is no harm in taking 

(2.8) 



__ T) 1 



where n > is a real number and A is a constant with mass dimension 
2. 

Using (2.8) in (2.7) ,it is obtained that 



\2 Ca n 
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In the late universe, a(t) is large, so this equation reduces to 



where 



and 



C 



1 

67 



Ca" 
- 1 



Da r 



a 



18n/3' 



a 



D = — 

67 



16ttG 3/5 
1 



cr 



(2.9a) can be re-written as 

•■ ^ x a 2 Ca n+l 

a + (2-2ra) — = — 

a 



167rG 9/5. 

Da" 



1 - 



which integrates to 

2 E 2C r a r ' 



/ay _ 

W " a 6 " 4n nA 



Da 



2n 



C J 18n/3' 



3a 



. (6 - 3n) C(6 - 2n) J 9n/3(6 - 4n) 
where .E is an integration constant having dimension (mass) 2 



(2.9a) 
(2.96) 
(2.9c) 

(2.10) 
, (2.11a) 
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(2.11a) is the modified Friedmann equation (MFE) giving cosmic 
dynamics. Terms, on r.h.s. of this equation, are imprints of curvature. 
The first term of these, proportional to a~( 6_4n ) emerge spontaneously. 
It is interesting to see that this term (the first term on r.h.s. of (15a)) 
corresponds to matter density if n — 3/4 i.e. for this value of n, it 
reduces to Ea~ 3 and yields the density of non-baryonic matter being 
spontaneously induced by curvature. So, we recognize it as dark matter 
density. 

Thus, for n = 3/4, (2.11a) looks like 



E 16a 



729/5 



+ 



32Ca 3 / 4 
A5A 



6C \' 



(2.116) 



On r.h.s. of (lib), there are terms proportional to a 3//4 and a 3 / 2 . If 
the density term pd e = 4Ca 3 / 4 / '15 An G is put in conservation equation 



p dc + 3iJ(p dc + p dc ) = 
with pde = wpdc, we obtain 



w — — - < -1. 
4 



(2.11c) 



(2.11d) 



This result shows that pa c = 4Ca 3 / 4 /IhAnG behaves as phantom dark 
energy density being induced by f(R)— gravity. 
Now, (2.11b) is re-written as 

2 



H 2 
with 



8nG 16a 
-3-Pdm + 



Pdn 



I 32Ca 3 / 4 


r 5Da 3 ^i 




I 1 - 6C I 


J + 45A 




3E 





' - 8nGa 3 ' 

Using current value of pd m as 0.23p° r , (2.11f) yields 

Pdm = 0. 



(2Ale) 



(2.11/) 



(2.12a) 
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where ao = a(to), 3E/8nG = 0.23p° r aQ and 

= 3H 

with H = lOOkm/Mpcsec = 2.32 x 10" 42 GeV/i being the current 
Hubble's rate of expansion and h = 0.68. The present age of the 
universe is estimated to be t ~ 13.7Gyr = 6.6 x lO^GeV" 1 [19]. So, 



Hq 1 = 0.96t - 



Further, ao is normalized as 



a = 1. 



Connecting (2.12a) and (2.12c), it is obtained 

_ 0.23/£ 

Pdm o 



(2.126) 



(2.12c) 



(2.12d) 



WMAP [19] gives decoupling of matter from radiation at red-shift 

1 



Z d = 1 = 1089. 

ad 



(2.13) 



So, it is supposed that dark matter begins to dominate cosmic dynamics 
when a > 

Now, (2. lie) is re- written as 



a) -^K Pdm+ 243^G/ +Pd n 1 - 2A/ 



where 



Pdc = PX /4 



with p° dc = 0.73p° cr = 4C/15AttG using a = 1 and 

3CK 



A = 



5D ' 



(2.14a) 



(2.146) 



(2.14c) 



where C and D are given by (2.9b) and (2.9c) respectively. 



ACCELERATION AND DECELERATION ... 11 
The Friedmann equation (FE) (2.14a) is obtained from f(R)— grav- 
ity. Comparing it with general relativity based FE 

d\ 2 8nG A 



it is obtained that the constant term in (2.14a) behaves as cosmological 
constant 



16a , n _ , IX 

A = (2.14d) 

243/5 V ; 



with vacuum energy density 



PA = = 243^G- (2 ' 14e) 

It is interesting to see that (2.14c) contains a term — (pdc) 2 /2A anal- 
ogous to brane-gravity correction to the Friedmann equation (FE) 
for negative brane-tension [13] and modifications in FE due to loop- 
quantum effects [14]. Here A is called cosmic tension [9, 10]. (2.14c) 
shows that cosmic tension A depends on coupling constants a and (5 
in the gravitational action (2.1). Moreover, positive cosmological con- 
stant, too, emerges from curvature. 

From (2.12a) and (2.14b), it is found that pdm ~ a~ 3 and pdc ~ o?^- 
So pdm decreases and pd c increases with expansion of the universe. So, 
it is natural to think for values of these to come closer and to be equal 
at a certain time t*. At this particular time, we have 

0.23a,- 3 = 0.73a 3 / 4 

using (2.12a), (2.12c) and (2.14b) as well as a* = a(t*). This equation 
yields 

/23\ 4 /is 

a* = (73) • (2.15) 
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It shows that , for a < a* = (23/73) 4/15 , p dm > p dc > p 2 dc . So, (2.14a) 
is approximated as 



a\ 2 



8ttG 



Pdm + 



2a 



243/3ttGJ 



\2 0.23# 2 


H 


ah 


I ~ a 3 




' ~B. 



Connecting (2.12a) and (2.16), it is obtained that 



with 

obtained using (2.14e). 
(2.17a) integrates to 



B 



p° dm 10A8(3H 2 



PA 



a 



(2.16) 



(2.17a) 



(2.176) 



a = B^sinh 2 ' 3 



l/3\3/2 



(2.18a) 



This result is not consistent with the scale factor, obtained in standard 
model of cosmology during matter dominance. So, to have a viable 
cosmology, it needs to be approximated as 

2/3 



a = a d 
which is possible till 

2.B 3 / 2 



i+ (t-u) 



2a 



3/2 



{t-t d )+sinh-\a d /Bf' 2 <l 



(2.186) 



(2.18c) 



as sinhl = 1.18 ~ I. Here a d is given by (2.13), which is the scale factor 
at time t = t d = 386kyr = 2.8 x 10'%. Connecting (2. 13), (2.17b), 
(2.18a) and (2.18c), it is obtained that 



a 



f3H 2 ~ 



> 10" 



(2.18d) 



Connecting (2.16e) and a/ f3Hl ~ 10 6 from (2.18d), it is evaluated 
that 

p A = 2.19 x 10~ 8 Pcr = 5.48 x 10~ 55 GeV 4 . (2.18e) 
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The approximated form (2.17a) is obtained when a « a*, but, as 
discussed above, pd m — Pdc when a < a*. So, in the narrow strip around 
a = a* for a < a*, (2.14a) is approximated as 

2 87rGrr„ 2a 



+ 



243/3ttg] 
using p dm ~ p dc > p^ in (2.14a). 

Connecting (2.12a) and (2.19), it is obtained that 

3i 



(2.19) 



a\2 0.46# 2 







1 + 



a v 
a*/ 



(2.20) 



with = a/167.67/9ifQ. It is because at a = a*, mater-dominated 
phase ends and dark energy dominance begins. 
(2.20) integrates to 



a 3/2 = al /2 sinh 
being approximated as 



3H V0A6 



2a 



3/2 



(t — t d ) + sinh l {a d /a if ) 



3/2 



a = a d 



1 + ^Wp (t _y l2! 



2a 



3/2 



like the above case. 

Using (2.21), it is obtained that 



This result yields 

t* ~ 0.59t 

using values of a* (from (2.15)), Hq 1 (from (2.12b)) and t d . 



(2.21) 



(2.22) 



It is discussed above that for a > a*,pd m < Pde- In this case, (2.14a) 
is approximated as 

'd\2 8ttC7 



(:)" 



-Pd, 



(2.23) 
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Thus, in the late universe, a phantom model is obtained from curva- 
ture without using any source of exotic matter. But this model contains 
a correction term — AirGp\ e /3\ due to curvature-induced cosmic tension 
A being evaluated below. 

Connecting (14b) and (2.23), it is obtained that 

_ 3/4 _ 073p° r 



(2.24) integrates to 
a(t) = 



0.73# V /2 



2A 



■}]■ 



(2.24) 



0.73 P ; 



2A 



+ 



1.26 - 



0-73pg r 
2A 

2-i -4/3 



--ff o V0J3(f-f.)} ] (2.25a) 

as a* 3 ^ 4 = 1.26. (25a) shows that phantom model, obtained here, is 
singularity- free. 

From (2.25a), it is obtained that 

-1.7 pi 11 



a = 0.27H$a 5/2 



A 



a 

3 



-3/4 



This shows a > , when 



Further, (2.25a) yields 



1 = a = 



A 3 ' 



(2.256) 



(2.25c) 



0-73p°. 
2A? h 



+ 



1.26 - 



0.73pO r 



2AP h 

3 > 1 2 l - 4 / 3 

--//oV0J3(to-i*)} J 



(2.26) 



(2.28) 



Using (2.22) for t* in (2.26), A is evaluated as 

A = 5.77p° r . 

(2.25a) exhibits accelerating universe when t > t*. Thus, a transition 
from decleration to acceleration takes place at 

t = U = 0.59t (2.29a) 
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and red-shift 



1 /73\ 4 / 15 



(2.296) 



which is within the range 0.33 < < 0.59 given by 16 Type supernova 
observations [2]. (2.25a) shows that universe expands till pde becomes 
equal to 2A as it grows with expansion. It happens till a{t) increases 
to a pc satisfying 



•-pe 



3/4 



2A 



0.73p c 

Thus, expansion (2.25a) stops at time 



1.26- 



0.73p c 
2A 



(2.30) 



(2.31) 



3. Cosmic energy conditions as well as acceleration and 
deceleration during phantom era 

In what follows, it is found that correction term, — 47rGp^ e /3A due 
to curvature- induced cosmic tension A in (23), change the bahviour of 
the phantom model drastically. 

From (2.11c) and (2.23), it is obtained that 

- = 5- 3(p + p) 1 - - 

a 6 1 LA 

The correction term, in this equation, is caused due to curvature- 
induced cosmic tension A in (2.23). This type of equation was obtained 
earlier in [16, 17] in the context of RS-II model in brane-gravity. 



2p{l 



_P_ 

2A 



}]■ 



(3.1) 



In the GR-based theory, (3.1) looks like 

a _4nG, 
a 



■[p + 3P]. 



(3.2) 



Comparing (2.5) and (2.6), the effective pressure density P is given 



as 



p + 3P = 3(p + p) 



1 - 



P 
AJ 



-2p 



1 - 



_P_ 

2AJ 



(3.3) 
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Using (2. lid), (3.3) yields the effective pressure density in the cur- 
vature induced phantom model as 



(3.4) yields 



This equation shows that WEC is violated when p < 3/7A = 2.47p° r 
with A given by (2.28). Moreover , p + P = for p = 2.47p° r and 
p + P > for p > 3/7A = 2.47p° r . 

From (2.14b), it is found that energy density p for phantom fluid 
increases with increasing scale factor a(t). It is interesting to note from 
(3.5) that curvature-induced phantom fluid, obtained here, behaves 
effectively as phantom violating WEC till p < 2.47p° r , but phantom 
characterisic to violate WEC is suppressed by cosmic tension when p 
increases and obeys the inequality p > 2.47p° r . 

Further, using (2. lid) in (3.3) , it is also obtained that 

P + 3P = -^X -4A- (3 ' 6) 
Connecting (2.28) and (3.3), it is obtained that SEC is violated when 
p < 11A/7 = 9.07p° r . Also, it is found that p + 3P = for p = 9.07p° r 
and p + 3P > for p > 9.07p° r . 

Thus, it is obtained that (i) WEC is violated for p < 2.47p° r , (ii) for 
2.47p° r < p < 9.07p° r WEC is not violated, but SEC is violated and 
(hi) for p > 9.07p° r neither of the two conditions is violated. Also it is 
interesting to note that these corrections cause effective phantom divide 
at 

p = 2.47p° r . (3.7) 
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Moreover, these results suggest that a transition from violation of 
SEC to non-violation of SEC will take place at 

p = 9.07p° r . (3.8) 

Also, universe will super-accelerate till p* < p < 2.47p[! r , accelerate 
when 2.47p° r < p < 9.07p° r and decelerate when 9.07p° r < p < 11.54p° r 
as expansion of phantom phase of the universe will stop at p = 11.54p° r . 

These results are also supported by (2.25b) as (2.25b) and (2.14b) 
yield 



a = 0.27# n 2 a 5/2 



0.29 



2.68p 



o 



P 



(3.9) 



Connecting (2.14b), (2.25a) , (2.28) and (2.29a), it is obtained that 
p = 0.73p° r [0.06 + {1.094 - 0.32H (t - O^o)} 2 ]" 1 . (3.10) 

(3.7) and (3.10) yield that effective phantom divide is obtained at 
time 

t~2.42t . (3.11) 

(3.8) and (3.10) yield that transition time for violation of SEC to 
non-violation of SEC will take place at 

t~3.44t . (3.12) 

These results imply super-acceleration during the time interval 0.59to < 
t < 2.42t , acceleration during the time interval 2.42t < t < 3.44t 
and deceleration during the time interval 3.44t < t < 3.87t - Expan- 
sion, driven by phantom, will stop at time t = 3.87to as pa c will acquire 
the value 2A by this time. 

When t > 3.87to, deceleration, driven by matter, will resume and 
Freidmann equation reduces to (2.17). 
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4. Cosmic collapse using classical approach and its quantum escape 
4(a). Re-appearance of matter-dominance and cosmic collapse 

As mentioned above, dark energy terms are switched off, in (2.14a), 
at p dc = 2X = 11.54p° r (obtained from (2.28)), a = a pc given by (2.30) 
when t = t pc = 3.87t given by (2.31). So, for t > t pc = 3.87t . (2.14a) 
will look like 

o^^m^-itn <««> 

where 

pZ = 3.68 x 10- 6 p° r , (4.16) 

Pa 2a K J 

which is evaluated using (2.18d), (2.30) and (4.1b). and 

as 

pll = 2A = 11.54p° r . (4.1e) 

(4.1a) integrates to 

a(t) = A 1/3 a pc sinh 2/3 sinh' 1 ^ 2 + 2.22 x \Q- A H Q {t - i pc )] (4.2a) 
yielding 

- = -q[2.22 x 10" 4 # ] 2 [cosec/i 2 #(t) - 2] (4.26) 

with 0(i) = sm/i _1 A- 1/2 + 2.22 x 10- 4 # (t - t pc ). (4.2b) shows decel- 
erated expansion cuased by matter dominance till 

sinh6(t) < 1/V2 = 0.707. (4.2c) 

The result (4.2a) is obtained when expansion is driven by the term 

3 r - . / a \ 3- 
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in FE. 

Moreover, though at a = a pe , 

"*(^) 3/4 { 1 -(^) 3/4 } 

vanishes, it will be negative for a > a pe . So for a > a pe , Friedmann 
equation (2.14a) is obtained as 



a \ 3 / 4 |Y a \3/4 

1-1 



(4.3) 



Obviously, the negative terms in (4.3) will try to stop expansion and, 
on sufficient growth of a(t) upto a m , expansion will reach it maximum 
such that a a =a m = and a m satisfies the equation 

■ f d m ^ 3/4 |- i d m ^ 3/4 

^pe ' 

being approximated as 



Pdm- 



(4.46) 



1 A-^2Af^V / {(^Y -l}. 

Va pc / lva pe / J 

(4.4b) yields the solution 

(^)-{5[ i+ ^^r- i+4 - 25xir7 < 44c > 

The negative sign (— ) is ignored here as it yields a m < a pe , which is 
not possible. Using (4.4c), it is obtained that 



Pdm 



pe 
Pdm 



1 - 1.28 x 10" 



(4.4d) 



Using (4.4c) in (4.2a), It is obtained that the time t — t m corre- 
sponding to a = a m is given by 

sinh^h.- 1 ' 2 + 2.22 x 10- 4 H (t m - t pc ) 

1 



sink 1 



l + \/l + 2p^/A 



-4-s/A 

sin/i _1 (0.18123) = 0.18123. 
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(4.4e) confirms deceleration during time period t pe < t < t m as it 
satisfies the conditon (4.2c). 

So, for t pc < t < t m , (4.2a) is obtained as 

~ 1 2/3 

a(t) = a pc 1 + 2.22 x 10- 4 # A 1/3 ap e 3/2 (* - *pe) • (4.5) 
Connecting (2.12b), (4.1c) and (4.4e), t m is evaluated as 

t m - tp C = 5.32 x 10" 4 t = 3.51 x 10 38 GeV _1 = 694.4kyr. (4.6) 



Further, it is interseting to note that the curve a = a(t) will be 
continuous at t — t m , but the direction of tangent to this curve (pointed 
at a = a m ) will change because it will attain its maximum at t — t m 
yielding d < for t > t m , which used to be positive for t < t m . It means 
that universe will retrace back at t — t m and will begin to contract. 
During the contraction phase, term proportional to a~ 3 will dominate 
over terms proportional to a 3 / 4 and (4.3) will yield 



1/2 

\a m ) /. 

On integrating (4.7), it is obtained that 

a(t) = A 1/3 a m {sinh[2.22 x 10' 4 H }(t col - t)} 2/3 , 

where 

icoi = t m + [2.22 x 10~ 4 Ho}" 1 sink' 1 A~ 1/2 
= t m + 3.33 x 10 2 t . 



(4.7) 



(4.8a) 



(4.8b) 



(4.8a) shows that a(t) = at t = t co \ and energy density of the 
universe 

3 / a \3/4i 



p = 

- Pdm 



J 



, a \ 3 



(4.9) 
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will diverge. It means that universe will collapse at t — t co i. 

This result, suggested by the classical mechanics, is unphysical due 
to divergence of cosmic energy density. So, to have a viable physics 
around the epoch t = t co \, the diverging component of the density 
of the universe, which is proportional to a~ 3 , needs to be finite. As 
this unphysical situation is predicted by classical mechanics, we have 
no other alternative other than to resort to quantum gravity. In what 
follows, we proceed in this direction, which is justified due to extremely 
high energy density and large curvature in the small interval of time 
near this particular epoch t co \. This situation is analogous to very 
early universe, where quantum gravity effects are dominant. Earlier 
also, quantum gravity corrections were made in the equations of future 
universe under such circumstances to avoid finite time singularities 
[15, 16, 17]. 

4(b). Further analysis of contracting universe through 
classical theory approach and conformal scalars 

(4.9) indicates that energy density of the contracting future universe 
will increase very high with increasing time. As a consequence, stellar 
equilibrium will break and, on sufficient rise in cosmic energy density, 
all matter in the universe, including celestial bodies, will get smashed 
to highly-relativistic baryonic as well as non-baryonic elementary par- 
ticles. Moreover, binding energy of celestial objects will be released 
as radiation with high temperature being in the high energy states 
(though , here, universe is driven by dark energy and dark matter as 
these constitute dominant components of the cosmic energy density 
and emerge from gravitational sector on which model is based, but 
presence of non-dominant baryonic comonent is not ignored). Energy 
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distribution of these elementary particles (bosons and fermions) is given 
by 

Pr = Pb + Pf = — \g b + T 4 (4.10) 

in natural units used in this paper. Here Pb(Pf) is the energy density 
and gb{gf) is the helicity for bosons (fermions). T is the temperature. 
Moreover, pressure density for these particles, is obtained as 

Pr =Pb+Pf = ^[Pb + Pf] = ^p. (4.11) 

It shows that, on growth of energy density in the universe upto suffi- 
ciently high level, energy distribution will also change causing a change 
in the EOS of dominating fluid from p = (for dark matter addressed 
in 4(a)) to EOS given by (4.11). Incorporating this change, cosmic 
energy density is obtained from conservation equation 

p r + 3H(p r + p r ) = 

as 

Pr = P^)\ (4-12) 

where p c h and a c h are energy density and scale factor at the epoch 
t = t c h (when highly relativistic particles and radiation are produced 
due to smash of cosmic matter at the epoch). 

EOS p r = 3p r for an ideal fluid is obtained when trace of the energy- 
momentum tensor components vanish. At the classicl level, this is true 
for conformal scalars <fi also. So, it is reasonable to think of existence of 
conformal scalars too alongwith other fields repesenting bosons and 
fermions giving energy density like (4.12). The action for conformal 
<f>(x, t) Is given as 

S^ = J d 4 xV=g~l [g^d^d^ - l -R<f\ . (4.13) 
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According to above arguments, for t > t c h 5 cosmic matter will break 
to elementary particles having different energy distribution (4.12), so 
energy distribution (4.9) will be replaced by (4.12) in FE as source of 
Pdm ~ a~ 3 will vanish. So, for t > t c h, (4.7) and (4.8a) will not be 
valid. Now, in this stage, (4.3) is modified as 

87TC7 



with 



A ch = Pch/PA- 



(4.14a) yields 



8ttC7 



-Pch 



i + A-(-iLY 



1/2 



(4.146) 



(4.14c) 



1/2 



replacing (4.7). 

On integrating (4.14c), it is obtained that 

a(t) = A 1 / 2 a ch {sinh[ ] f^^(C - t)] } 
where 

sinh^K' 1 ' 2 . 



t 



ch 
col 



(4.15a) 



(4.156) 



32irGp ch j 

Like (4.8a), (4.12) and (4.15a) also exhibits cosmic collapse in fu- 
ture as a(t) — > and p — > oo when t — > i^oi- Thus, in the changed 
circumstances too, classical theory predicts cosmic collapse in future. 

All these results indicate that contracting phase of the future uni- 
verse will undergo a revival of states of the early universe, where ele- 
mentary particles in thermal equillibrium will dominate. According to 
the standard model of cosmology, thermal equilibrium of these parti- 
cles are maintained at and above 10 _3 GeV. So, using it in (4.10), it is 
obtained that 

11tt 2 . 

(4.15c) 



p ch = H^! T 4 = Lg x 10 -12 GeV 4 



60 
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4(c). Quantum escape of collapse 

In what follows, a possibility to escape cosmic collapse, obtained by 
classical appproach, is probed using quantum corrections. 

Quantum field theory suggests that anamolous terms arise in energy- 
momentum tensor components of conformal scalar <fi [16, 17, 20] and its 
trace does not vanish contrary to the results obtained through the clas- 
sical approach ( as mentioned above, trace of the energy -momentum 
tensor vanishes for conformal scalars according to the classical theory) . 
These terms are advantageous being free from ultra-violet and infra-red 
divergences with respect to energy modes. So, quantum gravity yields 
anamolous terms as renormalized energy density pa and pressure den- 
sity pa as 

Pa = [80^p [ 3g4 + 6 ™ + 18i/2 ^ _ (4 ' 16a) 

and 

P * = -^< W4+8 ^)-l8o|bi 2S+12i///+18ff2//+9i/2 l- 

(4.166) 

with finite positive integers N being the number of conformal scalars. 
On making quantum corrections, (4.14a) is obtained as 

= (p r + p A ) + p\ (4.17) 

using (4.14b) and p\ as p\ after quantum corrections. 

Due to correction terms in (4.17), a(t) (being the solution of this 
equation) will be different from (4.8a) or (4.15a). But, the required 
solution is not possible unless pa is obtained as a suitable function of 
a(t). Under these circumstances, there is an alternative approach to 
assume a solution and find conditions to make this soluition capable 
to yield physically viable results and satisfy (4.17). (4.8a) and (4.15a) 
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show that diveregence is caused by dominating component of cosmic 
energy density when a — > as t — > Now, if quantum corrections 
remove this divergence, a(t) should be finite at t — t^. In cosmology, 
mostly a(t) is found either as a power-law solution or as an exponential 
solution of FE. According to investigations here, t = t^ ch is the epoch 
where this universe is expected to end. So, if a(t) ~ (t^i — t) n for the 
contracting universe for t > i c h, diveregence can not be removed. It 
can be checked from (4.16a) that pa itself will diverge for this form 
of a(t). So, these arguments suggest to choose an exponential form of 
a(t) during the quantum regime as it is possible for this form to yield 
finite value a co i = 0(^1)- 

With this motivation, for t > t c , solution of (4.17) is taken as 

a\t) = a col exp[\{D(C - t)}\ + C^D^ - f)} 7 ] (4.18) 

with 7 being a positive real number, C being an arbitrary dimensionless 
constant and D being a constant of mass dimension used for dimen- 
sional correction. Here t — t c is the epoch near t = t^ ch , when quan- 
tum gravity effects begin to dominate. In what follows, it is shown that 
(4.18) is physically valid satisfying (4.17) and numerical values Ci, D, 7 
and a co \ are determined. 

From (4.18), it is obtained that 

H = -D- 7 C7 1 D[D(**-0] (7 - 1) , 

H = 7(7-1)^^-0]^, (419abcd) 
H = - 1 ( 1 -l)( 1 -2)C 1 D 3 [D(tt l -t)]^\ 
H = 1 ( 1 -l)( 1 -2)( 1 -3)C 1 D 4 [D(C-t)]^. 
From (4.16a) and (4.16b), we find that terms containing H and H 
are dominant in pa and pa- So, minimum possible value for 7 is taken 
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as 

7 = 2, (4.20) 
because H and H do not vanish even at t — for this value of 7. But 

H = = H. (4.21a, b) 

Also, at t = t^oi 

H = —D and H = 2C X D 2 , (4.22a, 6) 

which is evaluated using 

M£>(f£. " (C " e])] (7 - 2) = limM - [C + e])] (7 - 2) = 1 (4-23) 
for 7 = 2. 

Now, it is important to check whether a q (t), given by (4.18), satisfies 
(4.17) with p A from (4.16a). Connecting (4.16a), (4.17), (4.18) and 
(4.19a)-(4.19d), it is obtained that 

^M 2 D 2 [1 + AC^DiC -t)} + 4C 1 2 {D(^ 1 - t)} 2 ] = Pdm [l - 

4|{^(C " t)}\ + 4(2 + C^DiC - t)} 2 ] +p A + - 4C X 2 ) 

+8C 1 { J D(C - f)} + 24C 2 {D(C - t)} 2 } (4.24) 

approximating upto second order only as — t)\ < 1 , because, 

according to above arguments, quantum gravity is effective near the 

epoch t = i£oi- Comparing constant terms and coefficients of {-D(^oi ~~ 

t)} as well as {D(tf^ - 1)} 2 in (4.24), it is obtained that 

3 ND 4 
—MpD 2 = fr | fc<+PA + _ { _ ?( i_40»), 

3 87V _D 4 

-W = -4 ft | t-a + _ Cl , 
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(4.25a,b,c) 



Incorporating (4.16a) and (4.25a,b,c) in (4.18), it is obtained that 

3 



Pr 



Pi 



8nG 



PA 



3 9 9 iVL> 4 , 
= —MlD 2 + , - (1 
8tt p 60(4tt) 2V 



4^ 2 ). 



(4.25d) 



yielding (p r + p^) finite at £ = t^i- 

Planck scale is the fundamental scale. It suggests the largest energy- 
mass scale as Planck mass M P . At this scale, energy density is obtained 
~ 10 75 GeV 4 . As, at t — t^, finite p r is obtained through quantum cor- 
rections. As, at t — t^ l5 we obtain finite p r due to quantum corrections, 
so 



J r\t=fi 



-Mi 



(4.26) 



where r > 1 is the real number. 
(4.25b) yields 



Ci 
2tt 



MpD 2 



3 - 



N 



60ttJ 



-Ap 



r\t=t° 



From (4.26) and (4.27), it is obtained that 



1 

27 



N 



60ttJ 



Ci = 



4 

r 



and 



D = M P 

as -D has mass dimension and C is dimensionless. 
Using (4.26), (4.28b) in (4.25c), we obtain 



(4.27) 



(4.28a) 



(4.286) 



3 

27 



N 

607T 



- 1 



c( + -d + - = o. 

r r 



(4.29) 
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(4.29) yields the real value of C, when iV < 60n. So, we can take 

N = 188. (4.30) 
Connecting (4.28a) and (4.30), C is obtained as 

ft = (4.31) 

r 

(4.29) and (4.31) yield 

r = 2tt. (4.32) 

and 

d = -1.995. (4.33) 

Connecting (4.25a), (4.26), (4.28b), (4.30), (4.32) and (4.33), we 
obtain 

pi = 2.56 x 10 75 GeV 4 . (4.34) 

Thus, it is found that vacuum energy density p\, given by (4.34), 
is extremely high and rises by 10 129 compared to p\ given by (2.18e) 
on making quantum correction. It means that quantum gravity feeds 
cosmological constant so well that it gets extremely fattened. Thus, it is 
interesting to see that, during quantum regime, cosmological constant 
is extremely heavy, but it is almost negligible in classical domain. 

As a q (t), given by (4.18), is valid during quantum dominance and 
it will have its maximum at t — t c , if quantum corrections begin to 
dominate at this epoch. So, from (4.19a), (4.20) (4.28b) and (4.33), we 
obtain 

H = -M P + 3.99Mp[M P (tf ol - t c )\ = (4.35) 

yielding 

t c = t col - 0.25t P (4.36) 
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with tp = Mp 1 being the Planck time. Thus, quantum gravity correc- 
tions will be effective for 

\M P (t col -t c )\ < 0.25. (4.37) 

(4.18), (4.20), (4.34b) and (4.36) yield 

a c = 1.13a col . (4.38) 

(2.30), (4. lb), (4.4d) and (4.15a) yield 

a ch = 1.468 x 10~ 12 . (4.39) 

From (4.1c), (4.15a), (4.15c), (4.36) and (4.39) , we obtain 

a c = 1.875 x 10" 34 . (4.40) 

Further, (4.38) and (4.40) imply 

a co i = 1.659 x 10" 34 . (4.41) 

Thus, solution of FE (4.17) with quabtum currection is obtained as 

a q {t) = lM59xl0' u exp[\{Mp(tf ol -t)}\-1.995{M P (t c ^-t)} 2 } (4.42) 

(4.41) yields H < for t < t col and H = M P [1 -3.99{M P (t-t^}] > 
for t > t co i. It shows contracting universe for 0.25tp < t < 
and expanding universe for < t < 0.26tp during quantum regime. 
It means that future universe will pass through the narrow quantum 
bridge from pre-collapse to post-collapse siate avoiding cosmic collapse 
suggested by classical mechanics. 



5. Summary 

Here, cosmolgy of the late and future universe is obtained from 
f(R)— gravity, obtained by adding higher-order curvature terms R 2 
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and R 3 to the Einstein-Hilbert linear in scalar curvature R. It is ex- 
plained in the first section on introduction that f(R)— gravity cos- 
mology is different from f(R)— dark energy models [6]. Moreover, 
problems of f(R) — dark energy models, pointed out in [7], do not ap- 
pear in f(R)— gravity cosmology [8, 9, 10, 11, 12] as well as in the 
present model. Here, curvature scalar contributes to both geometrical 
and physical components of the theory. Thus, it plays dual role as a 
geometrical as well as physical fields, which was obtained earlier in [21]. 

Here, it is found that, in the late universe, for the red-shift z < 1089, 
dark matter term emerge spontaneously and phantom dark energy 
emerge as imprints of linear and non-linear terms of curvature. It 
is found that, during 0.36 < z < 1089, dark matter dominated and 
universe expanded with deceleration as t 2 ^ 3 . A transition from decel- 
eration to acceleration took place at z — 0.36 and at time t = 0.59t 
(to being the present age of the universe). This transition is caused by 
dominance of curvature-induced phantom dark energy over curvature- 
induced dark matter. Dark energy gives anti-gravity effect and phan- 
tom dark energy exhibits this effect more strongly due to violation of 
WEC. So, phantom energy puts high jerk causing super-acceleration. 

f(R)— gravity inspired Friedmann equation, obtained here, contains 
two terms (i) 87rGpd m /3 ( with pd m being the dark matter density) 
and (i) (87rGpd e /3)[l — pd c /2A] ( with pd c being the phantom dark 
energy density). Here A = 5.777p° r is called the cosmic tension. It 
is interesting to see that Friedmann equation, obtained here, contains 
a correction term —4iiGp 2 c /3\ analogous to such term in Friedmann 
equations from RS-II model of brane-gravity [13] and loop quantum 
cosmolgy[14]. This correction is not effective in the present universe as 
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Pdc 2A as well as it is ineffective till pd e << 2A. But, as pd c will 
increase in future with the growing scale factor a(t), effect of this term 
will increase. It is found that, on sufficient growth of pd e , the effective 
EOS does not violate WEC (which characterizes phantom) , but violates 
SEC. On further increase in p de , even SEC is not violated. Thus, 
in the beginning of the phantom phase, universe will super-accelerate 
during the period 0.59to < t < 2.42to, it will accelerate during the 
period 2.42t < t < 3.44t and, during the period 3.44t < t < 3.87t , 
universe will decelerate even during the phantom phase. Phantom- 
dominance will end when p = 2A = 11.54p° r at time t = 3.87t - As a 
consequence, re-dominance of dark matter will begin giving decelerated 
expansion. But, as universe will expand, growth of pd e ~ a 3//4 will also 
continue giving pd c > 2A. It causes the term (8irGp<i e /3)[l — pd e /2A] 
to switch over from positive to negative. Growth of this negative term 
will try to slow down expansion more rapidly. As a result, universe will 
reach a state, when expansion will stop and scale factor will acquire its 
maximum value in finite time t m = 3.87to + 694.4kyr. When t > t m , 
universe will retrace back and contract. Results, obtained here, show 
that contraction will continue for sufficiently long period 333to and will 
collapse at time t co \ = 336.87t + 694.4kyr, where energy density of the 
universe will diverge and scale factor will vanish. These results are 
obtained using prescriptions of classical mechanics. In this paper, it 
is probed further whether quantum gravity corrections can save the 
universe from the menace of collapse. 

It is argued, here, that as in the very small time priod close to t — t co \, 
energy density is extremely high and curvature is very large, quantum 
gravity effects will be dominant. It is analogous to the state of very 
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early universe, where energy density is extremely high and curvature 
is very large. So, quantum gravity has a crucial role. Motivated by 
these arguments, quantum corrections are used in the future universe 
near time of collapse and it is found here that expected collapse of the 
future universe can be avoided. Moreover, results show that contracting 
universe, in pre-collapse era will expand with super-acceleration in the 
post-collapse era. 

Interestingly, energy density due to cosmo logical constant p\, during 
quantum regime, is found 10 129 times pa during regime of classical 
mechanics. It means that quantum gravity nourishes the cosmological 
constant well and makes it very strong, otherwise it is very weak under 
the rule of classical mechanics. So, when p\ is out of quantum domain, 
it falls down immediately by 129 order. 

The state of the universe near the epoch t = t co \ can be realized 
as revival of the state of early universe in future due to high energy 
density, large curvature and dominance of quantum gravity effects. If 
it is true, study of this state of the future universe can unveil many 
hidden knowledge of the early universe. For example, here, results of 
the cosmological constant near collapse time explains largeness of p\ 
in the very early universe and its extremely small value in the present 
universe. It is important to see that this result is obtained without any 
fine-tuning. 
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